STABILITY OF ENTROPY SOLUTIONS FOR LEVY MIXED 
HYPERBOLIC-PARABOLIC EQUATIONS 



KENNETH H. KARLSEN AND SULEYMAN ULUSOY 



Abstract. We analyze entropy solutions for a class of Levy mixed hyperbolic- 
parabolic equations containing a non-local (or fractional) diffusion operator 
originating from a pure jump Levy process. For these solutions we establish 
uniqueness (L^ contraction property) and continuous dependence results. 
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1. Introduction 

The subject of this paper is uniqueness and stability results for properly defined 
entropy solutions of mixed hyperbolic-parabolic quasilinear equations appended 
with a nonlocal (fractional) diffusion operator. These equations take the form 

(1.1) dtu + div/(u) = div(a(u)Vw) + C[u], 

where u — u it,x) is the unknown, {t,x) £ Qt (0,T) xR"^, d> 1, and T > 
is a fixed final time. The operator C is the generator of a symmetric positivity 
preserving pure jump Levy semigroup e*^ on L^{M.'^). 
Equation (1.1) is subject to initial data 

(1.2) u(0,x) =wo(a;) e (L^ni°°)(M'*). 
In (1.1), 

(1.3) /=(/i,...,/rf)eI^i'°°(R;R'^) 

is a given vector- valued flux function, a = (a-ij) > is a given symmetric matrix- 
valued diffusion function of the form 

(1.4) a = cr"(a'^)"', cr'^eK'^^^, 1 < K < d. 
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More precisely, the components of a are ay = J2k=i '^ik'^jk for «, j = 1, . . . , d. We 
assume that the matrix- valued function a" — (cr^j.) : M ^ R'*^-^ satisfies 

(1.5) (T° e iyi'°°(R;M'*''^^). 

Observe that we do not assume the matrix a(-) to be strictly positive definite, so the 
operator div(a(u)VM) may be strongly degenerate, and hence the phrase "mixed 
hyperbolic-parabolic" is justified. 

In terms of its singular integral representation, the nonlocal operator C in (1.1) 
takes the form 

(1.6) C[u]{t,x)= / [u{t, X + z) — u{t, x) — z ■ '^ul\z\<i] 7r((iz), 

where the singular Levy measure ■n{dz) is a positive, cr-finite Borel measure on 
R'' \ {0} satisfying 7r({0}) = 0, 7r((i(-z)) = -7r(dz), and 

(1.7) / (|z|'l|,|<i + |z|l|,|>i) ^(dz) <oo, 

JR<^\{0} ^ ^ 

where we note that z can be replaced by a certain regular jump function j{z) easily 
throughout the analysis. A typical example is provided by taking 

Az) = '^\z\<idz, ag(0,2). 

This example corresponds to the fractional Laplacian := — (— A)^ on M"^, which 
can also be defined in terms of the Fourier transform as 

This definition is employed in [2'')] to prove (1.6) in this case. 

Nonlocal operators like Aq are examples of a pscudodiffcrcntial operator V with 
a symbol a{Lo) > such that Vv{lo) — a{uj)v{uj). The function e~*°'^'^^ is positive 
definite, and thus, by the Levy-Khintchinc formula, it can be represented as 

a{uj) = lb ■ (jj + q[uj) + (l — e~*"'" — iz ■ oj l|2;|<i(2;)) 7r(dz), 

Js.''\{o} 

where 6 g E'^ represents the drift term, q{uj) = Ylt j=i Qij^i^^j is a positive definite 

quadratic function representing the pure diffusion part (g(w) = gives raise to 
the usual Laplacian — A), and the Levy measure Tr{dz) accounts for the jump (non- 
local) part. In our setting of £, cf. (1.6), we assume b = and q = 0, i.e, we are 
dealing with a pure jump operator. The key point is that any pseudodifferential 
operator V is the generator of a Levy process which is completely characterized 
in terms of the triplet {b,q,TT{dz)). For more details about the Levy-Khintchine 
formula and Levy processes in general, we refer to [11, 28, 29, 30, 43]. 

Integro-partial differential equations, also known as nonlocal, fractional or Levy 
partial differential equations, appear frequently in many different areas of research 
and find many applications in engineering and finance, including nonlinear acous- 
tics, statistical mechanics, biology, fluid flow, pricing of flnancial instruments, and 
portfolio optimization. Many authors have recently contributed to advancing the 
mathematical theory for quasilinear and fully nonlinear partial differential equations 
that are supplemented with a fractional diffusion operator arising as the generator 
of a Levy semigroup, addressing questions like existence, uniqueness, regularity, 
formation of singularities, and asymptotic behavior of solutions. 

For results with reference to fully nonlinear equations, such as the Hamilton- 
Jacobi-Bellman equation, and the (in this context relevant) theory of viscosity 
solutions, we refer to [2, 4, 5, G, 7, 16, 17, 18, 27, 31, 32, 42, 45, 44, 46, 47, 48], see 
also [9, 10, 23] for some concrete applications to finance. 



FRACTIONAL DEGENERATE PARABOLIC EQUATIONS 



3 



More recently, a number of authors [1, 3, 12, 13, 14, 15, 25, 33] have studied 
questions regarding existence, uniqueness, regularity, and temporal asymptotics for 
quasilinear equations, such as the fractal Burgers equation 

(1.8) dtu + 9,(mV2) = -i-dl,)^u, 

and more generally multi-dimensional fractional conservation laws 

(1.9) dfU + divf{u) = AaU, 

where the parameter a is assumed lie in the interval (0, 2). Of course, the excluded 
case a = 2 corresponds to the already fully understood viscous conservation law 
dtu + div/(u) — Au, solutions of which are always smooth in t > 0. Regarding the 
less studied case a G [1,2), it was recently proved in [24, 37] that solutions of the 
fractional Burgers equation (1.8) are also smooth in i > 0. In the case a < 1 for 
the fractional conservation law (1.9) the order of the diffusion part is lower than 
the first order hyperbolic part, so we do not expect any regularizing effect to take 
place. Indeed, for the fractional Burgers equation (1.8) with a < 1 it is proved in 
[3, 37] that solutions can develop discontinuities in finite time. Consequently, one 
should employ a notion of entropy solutions for fractional conservation laws (1.9), 
i.e., weak solutions satisfying an additional entropy condition, to ensure the global- 
in-time well-posedness. This is well-known for conservation laws dtu + divf{u) = 0, 
cf. Kruzkov [3S] , and the well-posedness theory of Kruzkov was recently extended 
to fractional conservation laws in [1]. 

In recent years the theory of Kruzkov [38] has been extended to quasilinear mixed 
hyperbolic-parabolic equations of the form 

(1.10) dtu + div/(u) = div(a(M)VM), 

where / and a satisfy (1.3) and (1.4)-(1.5), respectively. Since the diffusion matrix 
a(u) is not assumed to be strictly positive definite, (1.10) is strongly degenerate and 
will in general posses discontinuous solutions. In the isotropic case (with a(-) being 
a scalar function) the first general uniqueness result is due to Carrillo [1!)], who 
developed an original extension of Kruzkov's method of doubling variables to prove 
his result, cf. [34, 35, 39, 40] for some additional applications of his techniques. The 
anisotropic case (a(-) being a matrix-valued function) was first treated by Chen and 
Perthame [22], who developed a kinetic formulation and established the uniqueness 
result using regularization by convolution. An alternative proof of the result of 
Chen and Perthame, adapting the device of doubling variables, was developed in 
[8], cf. also [21, 20, 41] some other papers dealing with the anisotropic case. 

The main purpose of this paper is to extend the uniqueness and "continuous 
dependence on the nonlincaritics" results of [S, 21, 20, 41] to fractional degenerate 
parabolic equations of the form (1.1). We introduce the notion of entropy solutions 
and state the main results in Section 2 . Sections 3 (existence), 4 (uniqueness), and 
5 (continuous dependence on the nonlinearities and the Levy measure) are devoted 
to the proofs of the main results. 

2. Notion of solution and main results 
For i — 1 , . . . , d and k = 1, . . . , K, define 

Jo Jo 
for any ip £ C(R). Given any convex entropy function 77 : K. — > R, we define the 
corresponding entropy fluxes g = (gj) : M — *■ M'' and r ~ (vij) : M ^ K.'^^'^ by 

q'{z)=^'{z)fiz), r'{z)=^'{z)a{z). 

We refer to {"q, q, r) as an entropy-entropy flux triple. 
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We now introduce the entropy formulation of (1.1)- (1.2). 

Definition 2.1. An entropy solution of the initial value problem (1.1)-(1.2) is a 

measurable function u : Qt — * K satisfying the following conditions: 

(D.l) u e L"°iQT), u e L°°(0, T; L^[M.'^)), 

d 

(2.1) ^a,.aH€i2(gT), k = \,...,K, 

i=l 

and 

(2.2) // / {u{t, X + z) - u{t,x)f TT {dz)dxdt < +00. 
J Jqt Jr-^XIo} 

(D.2) Fork^l,...,K, 

d d 

(2.3) ^a,.Crfc'''(w) = V'(w)5]c>,.a(w), a.e. m Qt and m L\Qt), 



for any ^ G 

(D.3) For any entropy- entropy flux triple {ri,q,r), 



« « d d 

/ / {jl{u)dt(p + q^{u)^^^(p+ r,j{u)dl^^^(p^ dxdt 
(2.4) '■^'^ ''^=1 

+ // 'q{u)C[ip]dxdt + / ?7(uo)((5(0, a:)da; > n" + 771", 



JJQt Jm-^ 
for all non- negative ip <E C^{[0,T) x W^), where 

K d 

n 

" k=l i=l 



// r,"{u)j2i^d.^.QM)\{t.x)dxdt, 



™" ~ // / ri"{u]z){u{t,x + z) — u{t,x)) ip(t, x) '!T{dz) dx dt, 

J Jqt "'K''\{0} 

and 

r-l 



77"(u;z)=: / {1 - T)7f{{l - T)u{t,x) + Tu{t,x + z))dT. 
Jo 

If in addition, 

( \z\7r{dz) e L\R''\{Q}), 

(2.5) I or 

[ \z\^TT{dz) e L^W^XiO}) and u e L"" {0,T; BV{R'^)), 

then we can drop the fractional parabolic dissipation measure m" and replace 
(2.4) by the simpler condition 

« « d d 

/ / {viu)dt'p + Y qi{u)dx^ip + Y rrj{u)dl^^^Lpj dxdt 

(2.6) '^^ »J=i 

+ // ri{u)C[(f]dxdt + / 77(^0)1^(0, x) da; > n". 

We remark that the chain rule (2.3) is automatically fulfilled when a(-) is a 
scalar or a diagonal matrix, cf. Chen and Perthame [22], and in this case we can 
drop (D.2) from the definition. 
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Starting from the definition of £ (cf. calculations in the upcoming sections), we 
can replace the term 



ri{u)C[ip] dxdt — m", 

J JQt 

occurring in (2.4) by 

ri{u) [ip(t, X + z) — ip{t, x) — V(p • z] 7r(dz) dx dt, 



+ / T]'{u)[u{t,x + z) - u{t,x)]TT{dz)dxdt, Vre(0,l), 

J Jqt J\z\>r 

This formulation of the nonlocal term is directly related to the formulation used in 
[1] for fractional conservation laws. 

Our first result is the expected contraction property (and thus the uniqueness) 
of entropy solutions. 

Theorem 2.1. Suppose f and a satisfy (1.3) and (1.4)-(1.5), respectively, and 
that the Levy measure 7r(c?z) satisfies (1.7). Then there exists an entropy solution 
of (1.1)-(1.2). Let u,v be two entropy solutions of (1.1) with initial data u|t=o = 
Mo e (L^ ni°°)(M''), v\t=o = vo e (L^ n L°°)(R''). For a.e. t € (0,r), we have 

(2.7) / {u{t,x) ~ v{t,x))^ dx < {uo--vo)^dx. 



Consequently, if uo < vo a.e. in R'^ then u < v a.e. in Qt, so whenever uq = vo 
a.e. in M'', then u = v a.e. in Qt. 

Finally, if we replace (2.4) by (2.6), then the L^ contraction principle (2.7) 
continues to hold provided (2.5) is effective; it is sufficient that (say) only v belongs 
to L°°(0, T; BV(R'^)) in the case J \z\^ ■n{dz) < oo. 

This theorem generalizes to the "non-local diffusion" case the result of Chen and 
Perthame [22]. The proof follows that of Bendahmane and Karlsen [S]. 

Regarding the last part of Theorem 2.1, assuming vq S BV(R'^) it follows from 
(2.7) that V € L°°{0,T;BV{R'^)), as required. 

Our second result, which is a refinement of the previous theorem, reveals how 
the entropy solution u depends on the Levy measure Tiidz), and the nonlinear fluxes 
/, a (i.e., it is a "continuous dependence" estimate). 

Theorem 2.2. Suppose / and a satisfy (1.3) and (1.4)-(1.5), respectively, and 
that the Levy measure 7r{dz) satisfies (1.7). Let u G L°°{0,T; BV{R'^)) be the 
entropy solution of (1.1) with BV initial data uq E (L^ L°° n BV){W^) and 
with a Levy measure of the form T:{dz) — m{z) dz for some integrable function 
m : M''\ {0} ^ R+. 
Replace the data set 

(/, a,7r,Mo), a = cr°(cr)''^, Tr{dz) — m{z) dz 

by another data set 

(/, a, 7r(dz), Uq)? a = cr°(cr°)*'', TT{dz) = m{z) dz, 

where f,a°',TT,rh satisfy the same regularity conditions as f,a°',TT,m and moreover 
Vq G {L^ n L°°){M.'^). Denote the corresponding entropy solution by v, and assume 
that V G C{[0,T]; L^(R'^)). Suppose u and v take values in a closed interval / C R. 
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+ 6*2^^ II cr" - cr" II , 
Vi/i,oo(7).Rd) II ll-L (7,* 



For any t £ (0,T), 
(2.8) 

- v(i,-)llLi(M'') 
< ho-'follLi(Rci) +Clt f - f 

+ C3Vt. i [ \zf \m{z) - m{z)\ dz ] + Cit I \z\\m{z) - m{z)\ dz, 

where the constants Ci, z = 1, . . . ,4, depend on the L°°{Q,T; BV{M.'^)) norm of u. 
This theorem generaUzes results in [20, 21] to the "fractional case". 

3. Proof of Theorem 2.1 (existence) 

Although a detailed version of the existence of entropy solutions to (1.1) is pre- 
sented in [30] , to motivate the entropy condition and to present a brief sketch, let 
us consider the following accompanying problem containing a uniformly parabolic 
operator depending on a small parameter p > 0: 

(3.1) dfUp + div/(up) = diy{a{up)Vup) + C[up{t, •)] + pAup. 

It is standard to construct a smooth solution Up to (3.1), for each fixed p > 0. 
Indeed, it can be done using the Galerkin method and the compactness argument, 
see Chapter 5 in [2(i] and [37]. 

As usual, the game is to pass to the limit as p ^ and identify the entropy 
condition satisfied by the limit function u. We will be brief in establishing the 
following estimates, since most of them are similar to the ones in [22] and we will 
assume uq G W^^^ n n L°°(R'^), for general uq e L'^(R'^) one can follow the 
approximation procedure presented in [22]. 

The following estimates can be established for sufRciently regular initial data: 

II"pIIl°-(Qt) - ^' l"p(*' ■)lBV'(R'i) - ^' 

\\up{t2, •) - Up{ti, OIIlHK'') ~* ^' 1*2 - *i| ^ Oj uniformly in p. 
Hence there is a limit u such that, passing if necessary to a subsequence as p — > 0, 

(3.2) Up —>■ u a.e. in Qt and in U'{Qt) for any p G [l,oo). 

Next, we derive an energy estimate. To this end, fix a convex function rj and 
define q, r hy q' = rf f , r' = if a. Multiplying (3.1) by rj' yields 

d 

(3.3) dtiiiup) + divg(up) = ^ dfj^ijiup) + C[ri{up] + pArj^Up) - Up 

where Vp = v^p + v'^+ i/^ consists of three parts: 

(i) the entropy dissipation term 

i^l := pAr]{up) - p?]'{up)Aup = pii'{up) |Vwp|^ ; 

(ii) the parabolic dissipation term 

d K / d x"^ 

^«Vy("p) - 'n'{up)A\Y{a{up)^Up) = 77" (up) ^ ^-.Crfc("p) ; 

jj = l fe=l \i=l / 

(iii) the fractional parabolic dissipation term 

vl= j ri"{up;z){up{t,x + z) - Up{t,x)f T:(dz), 
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where r]"{up; z) = (1 — r)?7"((l — T)up{t, x) + TUp{t, x + z)) dr. 

In deriving (3.3), the "new" computation is the one showing that the commutator 

^Viup)] - ?]'{up)C[up] 

equals Vp, but this foUows easily from Taylor's formula with integral reminder; 

r]{b) - r]{a) = i f (a) {b ~ a) 

(3.4) / /-i \ , 

■ ' ' (1-T)r/"((1 -r)a + r6)drj (5-a)^ 

Specifying ri{z) = ^ in (3.3) gives 

r / E (e^-.c'^wI dxdt<c 

Jo k=l \^=l J 

and 

d d 

(3.5) E^-^C1(«p)^E^-.C'l(w) inL2(Q^). 

i=l i=l 

From this wc easily sec, as in [22], that (2.1) and (2.3) in Definition 2.1 hold. 

Regarding the non-local operator £, the same choice for -q reveals that (2.2) in 
Definition 2.1 holds. Now set 

n(dz) := (^z\^ l|3|<i + |z| l|2|>i) 7r(dz), 

and note that H{dz) is a bounded Radon measure. Introducing the short-hand 
notation 

Dp{t,x,z) = ' ' ' ' dn = U{dz)(g)dx(g)dt, 

\^\ l|.|<i + Vk|l|.i>i 

(2.2) translates into Dp being uniformly bounded in ^^((o, T) x R'^ x (R''\{0}); d/x). 
Consequently, we may assume that there is a limit function D such that 

Dp^D in L2((0,T) X R'^ X (R'^\ {0});d/i). 
Let us identify D. To this end, fix a smooth function ip in C^{Qt) and observe 
^(t, X) Mt,- + ^)^upit,x) ^^^^^ 
Qt JR-'XIO} |z| l|2|<i ^|z|l|3|>i 

^{t,x + z)-^{t,x) ^^(^^^)n(dz)d:.dt. 



\{o} k|l|z|<i + ^A^\Mz\>l 

u a.c. in 

u{t, x + z) — u{t, x) 



Now, using that Up - — *> u a.c. in Qt, we conclude that 



Dp - ' ' in L^((0,T) X R-* X (R-* \ {0}); d^)- 

1|2|<i + vl^|i|^i>i 

We are now in a position to pass to the distributional limit in (3.3) to recover 
the desired entropy condition satisfied by the limit u = limp^o Up- Note that to 
interpret (3.3) in the sense of distributions we use the formula 

(3.6) / C[^{x)](l){x)dx^ I ^{x)C[(j){x)]dx, 



which holds for all sufRciently regular (say, C^) functions $,0 : R'^ ^ R. This 
relation is easily obtained by a change of variables (f, x, z) i— > [t^x + z, —z) and an 
integration by parts in x. 
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We claim that the entropy condition satisfied by the limit u = limp^o takes 
the following form: for any convex entropy function r] and corresponding entropy 
fluxes q, r defined by q' ~ rj' f ,r' = rj'a, 

(3.7) dtr,{u) + divq{u) < J2 d.,x,nj{u) + C[i^{u)] - n"^" - jti"^" 

i-J 

in the sense of distributions, where 



k=l \i=l I 



is the parabolic dissipation measure with respect to u and 

m^'^ = / ri"{u; z) {u{t,x + z) ^ u{t,x))'^ 7r{dz), 



'R<*\{0} 

is the fractional parabolic dissipation measure with respect to u. 
In view of (3.2), to verify (3.7) we only need to argue that 

liminf // Updxdt> [f {n"^'' + m"^'') dxdt. 
JJqt JJqt 

First, JJq^ v], dxdt > for each p > 0. Second, thanks to the weak convergence 
(3.5) and a standard weak lower semi-continuity result for quadratic functionals, 

liminf r f 77"(«p)E fE^-.C*'^-("p)l ^'^^'^^ 

>r [ '/'HE (e^-.c^5c(«) 1 vdxdt, 

•^0 -^R" fci \i=i / 
for all test functions ip E . Similarly, 

liminf // / rj" {up-, z) {up{t,x + z) — Up{t,x))'^ (pTT{dz) dx dt 
> / / / ri"{u; z) {u{t,x + z) ~ u{t,x))'^ LpTT{dz) dx dt, 



IQt "'R''\{0} 

for all test functions ip E . Combining, we deduce that (2.4) in Definition 2.1 
holds. This completes the proof. 

4. Proof of Theorem 2.1 (uniqueness) 
We shall need approximations 77^ (z) of the functions 

?7±(z) := (z)± =max(±(z),0), zeR. 
We build these by picking nondccreasing approximations sgn^(z) of 

+ , , Jo, if z<0, , J-l, if z<0, 

sgn + (z) sgn [z)-=<. -n ^ n 

I 1, if 2; > 0, 10, if z > 0, 

and defining 

?7±(z) ~ f sgn±(e)d^, ZER. 
Jo 

For example, we can take 

0, if z < 0, f-l, if z < -e, 
sgn + (z) = ^ sin(5jz), if < z < e, sgn^(z) = < sin(^z), if -e < z < 0, 

1, ifz>e. 0, ifz>0. 
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The functions -qf are and convex. Moreover, 

77±(z)5=S7/±(z), zeM. 

Observe that (77|'^(- — c))^^]^ is a family of entropies. Given these entropies, we 
introduce the corresponding entropy fluxes 

qf{z,c) = (77±)'(^-c)/'(CK, z,ceR, 



Clearly, as e ^ 0, 

g±(z,c)^g±(2,c) :=sgn±(z-c)(/(z)-/(c)), z,ceM, 
rf{z,c) r'^iz,c) := sgn='=(u - c){A{u) - A{c)), z,c G M, 

r 

where the (matrix-valued) function A{-) is defined by A{z) = / a(^) d^. 

^0 

Observe that {r]^{- — c),q^{-,c),r^{-,c))^^^ is a family of entropy-entropy flux 
triples, so choosing t] = rj^ in (2.4) yields 

d d 
rif{u-c)dt(p + ^qfiiu,c)d^^ip+ ^ rf^^^{u,c)dl^^^ip^ dxdt 

i—l — l 



+ 11 r]f (u — c) C[lp\ dx dt + I r/ {uQ — c)(p{0,x)dx 
> // (r?±)"(u-c)f^(^a,,a(^.))%dzd^ 



(t^e )"(u — c; z) (^(t, a; -I- z) — u(t, x)) tp n{dz) dx dt. 

J JQt "'R''\{0} 

Moreover, 

{vt)"{u -c;z)= ^ (1 - r)(77±)" ((1 - r)u(t, a;) + Tu{t, a; + z), c) dr 

(1 - r)(sgn^)'((l - T){u{t,x) ~c) + T{u{t, x + z) - c)) dr. 





To proceed, the following simple observations will be useful: 

• sgn J (z — c) ^ — sgn + (c — z) and r]~{z — c) ~ rjf {c — z); 

• q~{z,c) = 9+(c, z) and r^{z,c) = r+{c,z); 
. (r;7)"(z-c) = (r;+)"(c~z). 

Employing these observations, we can rewrite the "— " part of (4.1) as 

d d 

^^r^+{c-u)dtf + ^qt^i{c,u)d.^,ip+ ^ r+^^{c,u)^l^^^^'j dxdt 
1=1 jj=i 

+ // 77^(c — u)£[(p] -I- / T]^ {c ~ UQ)(f{0, x) dx 
^4 2) •^■^'^^ 

> [[ ivt)"ic-u)J2(i29-^^^kiu))%dxdt 



Q-i 



fe=l i=l 



(r]t)"{c — u; z) (it(i, x -|- z) — u(i, x)) (/9 Tr{dz) dx dt. 
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To establish the contraction property (2.7) we shall employ the doubling- 
of-variables device of Kruzkov [38]. Let u = u(t,x), v = v{s,y) be two entropy 
solutions as stated in Theorem 2.1. Moreover, let ip = ip(t, x, s, y) be a test function 
in the doubled variables {t, x, s, y). To simplify the presentation, we introduce the 
following notation (with V x+y being short-hand for V^; + Vy) 



Cx[(p] := / [(p{t, X + z,s,y) - ip ~ z ■ Vx(pi\z\<i] 7r(dz), 

Js.'^\{o} 



Cy[(p] 



'\{0} 



[(p{t, x,s,y + z) - (p - z ■ Vj;<^l|^|<i] TT{dz), 

(p{t, X + Z,S,y + Z) - ip - Z ■ Vx+yipllzlKl TT{dz), 



In the "+" part of (4.1) written the entropy solution u{t, x) we choose c = v{s, y) 
and integrate the result over {s,y), obtaining 



(4.3) 



7]+{u-- v)dtip + ^q^i{u,v)dx,ip+ r+^J{u,c)^l^^^ipjdxdtdyds 

2—1 — 1 

rif{u~v)Cx['p] dxdtdy ds + jjj i]^ {uq — v)tp{0, x, s,y) dx dy ds 

K d 



> 



1 1 {rj+y'{u - v) J2 (E (^)) ^^dx dt dy ds 



k=l i=l 



'\{0} 



{ri^y'{u{t, •) — V] z) {u{t, X + z) — u{t, x)) ipTi{dz) dxdtdy ds. 



Similarly, in (4.2) written for the entropy solution v{s,y) we choose c = u(t,x) 
and integrate over (i, a;), thereby obtaining 



(4.4) 



> 



rj+{u-v)dsip + Yq^^,iu,v)dy^ip+ ^ dxdtdyds 

i = l i,j = l 

i]f {u ~ v)Cy[(p] dx dt dy ds + JJJ ri^{u — vo)Lp{t,x,0,y) dx dtdy 

K d 



{4)"{u - i') E (E %Crfc(«)) ^ dx dt dy ds 



fe=l i=l 



{i]ty'{u — v{s, •); z) (v{s, y + z) — v{s, y)) ip Tr{dz) dx dt dy ds. 



Adding (4.3) and (4.4) yields 



(4.5) /timo(£) + -?"conv(£) + ^diff (e) + -^tdiff (e) + Anit(e) > -?"diss(e) + -ffdiss(£), 
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where 

-?timc(£) 
-^conv (^) 

/diff(£) 
-ffdiff (e) 



77^ (w — v){dt + ds)(p dx dt dy ds 

d 

(m, u) {dxi + dy^)ip dx dt dy ds 

i=l 
d 

51 + dl^y,)'P dx dt dy ds 

rif{u — v)(^Cx[y^] + ^y[y^]j dxdtdyds 



-^init(e) = /// r]^{uo-v)ip{0,x,s,y)dxdyds 



ilf{u — vo)ip{t, X, 0, y) dx dt dy 



^diss(£) 



ivtriu-v) 

K 



-ffdiss(£) 



E 



R''\{0} 



Lp dx dt dy ds 



ivT)"{u{t, •) - v; z) {u{t, x + z)~ u{t, x)y 



+ ir]t)"{u,v{s, ■);z) {v{s,y + z) - v{s,y)) 
X ifnidz) dx dtdy ds. 
In view of the incquaUty "a'^ > 2ab'\ we have /diss(e) > -^diss(e), with 

UUe) = 2 /W ivtn^-v)Y.Y. d.:,XUu)dy^Qk{v)v dxdtdyds. 

k=l ij = l 

Arguing exactly as in [cS], it follows that 
lim(^/diff(£) - 7diss(e)) 

(4.6) i-rfi- d 

- r+M{dl,^+2dly^+dl^^)^ dxdtdyds. 

Fix a small number k > 0, and let us split C into two parts 

£[(P] = [ [(Pit, x + z)- (Pit, x)-z- V(/)l|,|<i] 7r(dz) 

J\Z\<K 

+ / [(f>it, x + z)- (pit, x) ~ z- V01|^|<i] 7r(dz) 

J\z\>K 

=: £40] +£«[(/>], V0eC2, 

and similarly 

The corresponding splitting of /fdiff(£) is written 

^fdifl(£) = Iiditf,Ki£) + -^fdiff (s)- 
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We also need to introduce the operator defined by writing 

\J\z\». ) 

with similar definitions for £,y, and C^j^y. Observe that (3.6) continues to 
hold for all these operators. The function obtained by replacing CJ^ with in the 
definition of /fjjfj(e) will be named I^^^^ie). 
Clearly, in view of (1.7), 

(4.7) l^fdiff.K(£)| <C||i^VL,.,o.xo.i / |z|'4rfz)^"0, 



for some constant C independent of k and e. 
Let us analyze /f'^jfj (e). By (3.6), 

^diff(e) = WW (i-l [^^(« - «)] + i-y - v)^^dtdxdyds. 

Specifying a ~ u{t, x) — v(s, y) and b ~ u{t, x + z) — v{s, y) in (3.4) yields 
f^j{u{t,x + z) -~v{s,y)) ~i]+(u(t,x) -v{s,y)) 
(4.8) = (?y+)'(w(t, x) - v{s, y)) {u{t, x + z)- u{t, x)) 



+ {rj+y'{u{t, •) - v; z) (u{t, x + z)- u{t, x)y . 

Similarly, taking a ~ u{t, x) — v{s, y), b ~ u{t, x) — v{s, y + z) in (3.4) yields 

■r]t{u{t, x) - v{s, y + z))~ 77+(u(t, x) - vis, y)) 

(4.9) = -{vtYi'^i'ti ^) - '"{s, y)) {v{s, V + z)- v{s, y)) 

+ {vt)"{u - v{s, ■);z) {v{s,y + z) - v{s,y)f . 

Adding the first term on the right-hand side of (4.8) to the first term on the 
right-hand side of (4.9) yields 

(?/+)'(ii(t, x) - v{s, y)) {u{t, x + z)~ u{t, x)) 

- (vtYi'^i'ti ^) - (s> y)) {v{s, y + z)- v{s, y)) 

= {''lt)'{'^{t, x) - v{s, y)) (u(t, x + z)~ v{s, y + z))- {u{t, x) - v(s, y)) 

< r]^{u{t, x + z)- v{s, y + z))- r]^{u{t, x) - v{s, y)), 

where we have used the convexity of 1]^ to derive the last inequality. 
In view of these findings, we can rewrite /f'Jiiff(e) as follows: 



/f«,iff(£) - ImJe) < //// [r?+(u(t, •) - i^s, •))] ifdtdxdyds 



where 

^fdiss(£) 



\z\>K 



(77+ )"(w(i, •) - v; z) {u{t, x + z) - u{t, x)) 



+ {-qtyiu ~ v{s, ■);z) {v{s,y + z) - v(s,y))'^ 



X (fni^dz) dx dt dy ds. 
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Consequently, 

-ffdiff (e) - ^fdiss(e) < JJJI Vti^ - ^O-CS+J^] dt dx dy ds, 

The next step is to first send k and then e ^ 0. Related to this, observe 
that 

lim -^fdiff (e) Ifdis (e) , Hm /fdi.s (e) = /fdiss (e) 

for each fixed e > 0, by the dominated convergence theorem. Moreover, we clearly 
have Htji^^o ^x+y[^] — ^x+y[^]- In view of this and (4.7), we conclude that 

(4.10) /fdiff (£) - -ffdiss(e) < JJJJ Tj^iu - v)C[ip] dtdxdyds. 

By (4.6) and (4.10), It follows from (4.5) and sending e ^ that 
(4.11) 

d 



{u - v)+{dt + ds)(p + ^q+(u, v)(a2;, + dy^)(p 



1=1 

d 



■ J2 ^tjiu^^)idl.,+2dly,+dly^)^ + ri+{u^v)Cx+y[^] dxdtdyds 



(uq — u) iy9(0, X, s,y) dx dy ds + JJJ ~ vq) ip{t, x, 0, y) dxdtdy > 0. 

Let us specify the test function ip = ip(t,x, s,y). To this end, fix a nonnegative 
test function (p = (f>{t,x) e Cj?°([0,oo) x R'^), and pick two sequences {6u}i,^o 
C^(0, i^), {<5^}^>Q C C^{B{0, ^)) of approximate delta functions, where B{0,fi) 
denotes the open ball centered at the origin with radius fi. Then take 

(4.12) ifiit, X, s, y) =9^{s- t)5^{y - x)(j){t, x). 

Simple calculations reveal that 

{dt + ds)(p = 9^{s - t)Sf,{y - x)dt(j){t,x), 
{dxi + dy^)(p = 9^{s - t)6f,{y - x)dxi(l){t,x), 
{dlx, + 2dly, + dly.)^ = e,is ~ t)S^iy ~ x)dl,^q^{t, x) 

and 

ip{t, X + z,s,y + z) - (p{t, X, s, y) 

= 0y{s ~ t)5^{y - x) {(j){t, x + z)- (i){t, x)) . 

Note that 0,^ = on (— oo,0] and so ifi{t,x,0,y) = 0. By the choice of the test 
function tp and the observations above, we deduce from (4.11) that 

JJJJ {'^ ^ v)'^9y{s — i)5^{y — x)dt4>{t, x) dxdtdy ds 

. d 

qf{u,v)0^{s - t)8^j_{y - x)dxi4'{t, x) dxdtdyds 



(4.13) 



d 

rfj{u, v)9u{s — t)S^{y — x)d'^.^. (j)(t, x) dx dt dy ds 

1,3 = 1 

[u — v)~^d,y{s — t)Sf^{y — x)C[(j)] dx dt dy ds + Iua,v{'^, m) > 0, 
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where 

/„o,.„(i/,^) := /// {uo- v)^0i,{s)6^{y ~ x)(l){0,x)dxdyds 

{uq - v)+ds {4>i,{s)5^,{y - x)(l){Q,x)^ dxdyds, 
with 

Ms)-^ OAr)dT^ 0,(r)dT'^l. 

Specifying ip ~ (j)^{s)5 ^j_{y — x) in the entropy inequality for v and noting that 9v{s) 
vanishes for s > i^, we obtain 



(uo - v)+dsip[s, X, y) dy ds 
(4.14) < //(mo - v)+6,{s)6^,{y - x)<^(0, x) dy ds + o{v) 

(mo -v)+5f,{y- x)(j}{{),x)dy ds, 



where the "o(i^)" term follows from an intcgrability argument. 
Hence, sending ly, i^i —^ 0, we deduce 

limsuphm sup/uQ^„(i/, fi) 



(4.15) < limsup 1 1 {uq- uo)+(5^(?/ - x)4>{0, x) dx dy 



{uo - vo) 4>{0,x) dx, 

with Uo = uo{x) and vq = vo{x). 

Keeping in mind (4.15) when sending fi, ly in (4.13), we conclude that 

JJQt\ 

(4-16) + ^ r+iu,v)dl,^c^+iu~v)+C[cj^]\ dxdt 

ij = l / 

+ / {uq- vo)+(p{0,x)dx > 0, 

where all the involved functions depend on {t,x). It now only takes a standard 
argument to conclude from (4.16) that Theorem 2.1 holds. Indeed, one chooses a 
sequence of functions < < 1 from C^{[0,T) x R'') that converges to l[o,t)xRd 
for a Lebesgue point t of Jjg^d(u — v)^ dx and then use the intcgrability of u,v to 
conclude the proof. 

Finally, let us prove the last part of Theorem 2.1, that is, we shall establish the 
contraction property for entropy solutions satisfying the simpler entropy condition 
(2.6) in which the fractional parabolic dissipation measure has been dropped. To 
this end, let us assume that (2.5) holds. Arguing exactly as before we arrive at 
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(4.5) without the /fdiss(£) term: 

' d 

{u - v)+{dt + ds)ip + ^ v){d^. + dy,)ip 



(4.17) 

+ (u — v)^{Cx['p] + ^yi'p]) dxdtdy ds 



+ {uq - v)^ (p\t=o dx dy ds + i i i {u - vq)^ ip\s=Q dx dtdy > Q 



The new treatment concerns the fractional term in (4.17) only, which we denote 
by J{y,iJL), i.e., J{v,^) := ff{u — v)~^{Cx['p] + Cy[ip]) dx dt dy ds; we employ the 
same test function (p as before, cf. (4.12). By letting z '—^ —z in the Cy term, 
keeping in mind that Tr(d{—z)) = —TT{dz), we obtain 

(4.18) 



J(j.,m)=//// / iuit,x)-vis,y))+e,is~t) 
Js.''\{o} 

X iSf,{y-x- z){(t){t, x-V z)- <f>{t, x)) 



Sf^{y — x)'S/(j>{t, x) • 2: l|2|<x Tr{dz) dx dt dy ds 



Sending the "temporal parameter" ly to zero in (4.18) yields 

{u{t, x)-v{s,y))+e,{s-t) 



J{n) := lini J(z/,^) 



'\{0} 

X iSf_,{y - X - z){(j){t, x + z)- (t){t, x)) 



~ ^ 2:)V0(i, x) ■ z l|2|<i J Tr{dz) dx dt dy ds 

Under the additional requirements listed in (2.5) we are also allowed to send the 
"spatial parameter" fi to zero, resulting in 

J := lim J(/^) 



(4.19) J J Jm^\{q} 



{u{t, x) — v{t, X + z))"*" ((/)(t, X + z) ~ (j){t, x)) 



— {u{t, x) — v{t, x))'^V(j){t, x) ■ z l\z\<_ij Tr{dz) dx dt. 

We divide the remaining discussion into two cases. 

Case 1: |2;|7r((iz) £ L^(M'' \ {0}). Adding and subtracting identical terms we 
obtain 

(4.20) J= hu{t,x) -v{t,x))+C[4>](t,x)dtdx + E, 
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where 



\E\ < I I / \v{t,x + z)-v{t,x)\\4>{t,x + z)-4>{t,x)\7T{dz)dtdx 



< 



so that we can employ the dominated convergence theorem to send 0^1, and 
consequently \E\ — > 0. 

Case 2: \z\^T:{dz) e L^R'^ \ {0}) and v e L°°(0, T; BV"(M'^)). In this case the 
error term \E\ in (4.20) can be estimated as follows: 

\E\ < (t ||«ILo.(o,T;W(K-) U"\\L^(Qr)) 1^1' e \ {0}; <dz)), 

and again we can employ the dominated convergence theorem. 
This concludes the proof of Theorem 2.1. 

5. Proof of Theorem 2.2 (continuous dependence) 

We again employ the doubling of variables device as in the previous section, but 
with a slightly different choice of the entropy function. For each e > 0, define 

-1, ifC<-e 
sgnaO = <( sin(^e), if ICI < e 
if ^ > £, 

which is a approximation of sgn (•). This choice gives rise to a approximation 
?7e(z) = agni;{^) d^ of the entropy flux \z\. As before, we introduce the corre- 
sponding entropy flux functions (u, c) , qf{u,c), and rf^ (u, c) . We now employ the 
doubling variables technique using the test function 

ip{t,x,s,y) e^{s - t)S^{y - x)Qa{t), 

where O^,, 5^ are symmetric approximate delta functions with support in 

and B{0,fj,), respectively. Fix a time t from (0, T). For any a > with < a < 

min(To, T — t), we define 

e„(i)-i/„(t)-i?„(t-r), H^it)^ f e,,{a)dG. 

so that e^(t) e^it) -9a{t~T). 

Proceeding as in the previous section (cf. also [2il]) and sending e — > 0, we find 

- JJJJ Kt - v\ 9t,{s - t)6^{y - x)Q'^{t) dxdtdyds < /conv - ^diff + -^fdiff, 
where 

^onv := JJJJ [G{u,v) - F{u,v)] ■ VxSfj.iy - x)9^{s - t)Qait) dxdtdyds, 

F{u, v) := sgn {u - v) {f{u) - f{v)) , G{u, v) sgn {u - v) {g{u) - g{v)) , 
r r r r ^ 

/diff:= //// Qa{t)9,{s~t)dl.,/^{y-x) I sgn{^ ~ v) 6%-' (0 dx dt dy ds, 

i.j = l 

K 

C(^) E {^tio^uo - 2<,(e)a^\(o + Ais.y.M ■ 

k=l 
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and Jfdiff = /fdiffi + /fdiffa with 

/fdiffi //// / \u^v\ e,{s - t)e„(t) 

JJJJ J\z\<l 

5^{y - X - z) - 5^{y - x) - V5i_,{y - x) ■ z 
X (m(z) — to(z)) (iz da: dy ds 



and 



-^fdiffa //// / |w - u|6'^(s - t)9Q(t) - a; - z) - - x) 
JJJJ J|z|>i L 

X (m(z) — m(z)) dz da; dt dy ds, 

By triangle inequality 

\u(t, x) — v{s, y)\ 6i,{s — t)5^{y — x)Q'^{t) dxdtdyds 

^ - nil \u{t,y)-v{t,y)\e,{s-t)Sf,{y-x)\e'Jt)\ dxdtdyds 
\v{t,y) - v{s,y)\ e,{s - t)Sf,{y - x) \e'^{t)\ dxdtdyds 
\u{t,x) - u{t,y)\ 6y{s - t)6^,{y - x) \Q'a{t)\ dxdtdyds 



— : £ + Rt + Rx- 

Keeping in mind that v G C{L^) and u £ L°°{BV), it is standard to show that 
lim Rt = 0, limsup \ Rx\ < Cfi 

and moreover, since also u{t) mq, v{t) vq as t ^ 0, 

limL ^ ||w(r, •) - OIIliir'') ^ \Wo - vo\\ l^r") ■ 
Following [2(j], using u E L°°{BV) we conclude that 

lim lim |/co„v| < Ct 11/ - gilLip(/) , 
and, exploiting also that / < C/fi, 

hm lim |/did < ||K - -^)(-" " -Th^^,.,^..., ■ 

It remains to estimate |/fdiff|. First, we nconsidcr /fdiffi- Using the Taylor and 
Fubini theorems we obtain 



|/tdiffj= /// / /' 

JJJ J\z\<lJO 



(1 - t)0^{s - t)Q^{t){m{z) - m{z)) 



\z\<l Jo 

X I / x) — w(s, y)! I?^Jp(y — .T — Tz) z • z (i.T j drdzdydsdt. 

Thanks to |u(<, •) — ?j(s, y)\ E BV{M.'^), an integration by parts yields 
(5.1) 

/fdiffi = III I I i^- - t)e^mm{z) - m(z)) 

JJJ J\z\<lJO 

X / W5f^{y — X — Tz) ■ z Dx {\u{t, x) — v(s,y)\) ■ z dx \ drdzdydsdt, 
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where the inner integral is taken with respect to the bounded Borel measure 
D {\u{t, ■) — v{s,y)\) ■ z. Since \D{u{t, ■) — v{s,y))\ < \D{u{t,-))\, the term inside 
the parentheses in (5.1), is upper bounded by 

/ / |V<5^(y -x~ Tz)\ \dD{u{t, ■))ix)\dy < jzp \u{t, ■)\bv(w^) I|V5,J|^wr.) , 

where we have used that \Du{t, •) | is finite and the Fubini's theorem to first integrate 
with respect to y. Hence, 



C f 
lim hm l/fdiffj < —t 



z|^ \m{z) — m(z)| dz, 

\z\<\ 



where C > is a finite constant. 

Similarly, relying again on the L°°(_B1/) regularity of w, it is not difficult to 
deduce via an integration by parts the estimate 

lim lim l/fdiffal < Ct \ \z\ |m(2;) — m{z)\ dz. 

ly^Oa^O J\z\>l 

Finally, we collect the bounds we have obtained so far and then optimize over fjL 
to obtain the desired continuous dependence estimate (2.8). 
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